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Section  I 

Summary  of  Accomplishments 
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I.  Summary  of  Accomplishments 

Under  this  Task,  the  following  problem  areas  have  been  addressed: 

1.  Gaussian  beam  (GB)  method  for  underwater  acoustic  propagation,  with 
emphasis  on  critical  reflection  and  head  wave  modeling  due  to  a  high 
velocity  bottom. 

The  objective  here  is  to  improve  the  narrow-waisted  paraxial  Gaussian 
beam  (PA-GB)  algorithm,  which  is  computationally  efficient  but  does  not 
account  for  the  total  reflection-head  wave  phenomena.  A  test  problem 
involving  a  line  source  in  a  homogeneous  ocean  above  a  homogeneous, 
semi-infinite,  fast  fluid  bottom  has  been  analyzed  in  detail.  Various  options 
for  augmenting  the  deficient  PA-GB  stack  have  been  proposed,  examined 
and  compared  with  a  rigorous  numerical  reference  solution. 

Details  of  the  analysis  and  results  obtained  so  far  are  given  in  Section  II. 
This  Investigation  will  continue  under  a  separate  Task. 


2.  Hybrid  ray-mode-PE  method  for  high-frequency  propagation  in  a 
tropospheric  surface  duct. 

The  objective  in  this  investigation  is  to  Improve  the  ray  prediction  code 
now  In  use  at  NOSC  by  self-consistent  hybrid  ray-mode-PE  parametrization 
In  those  observational  regions  where  ray  theory  fails.  The  canonical  model 
Is  a  line-source-excited  duct  with  laterally  homogeneous  bilinear  height 
profile.  An  exact  numerical  reference  solution  can  be  constructed  by  modal 
summation  for  comparison  with  the  results  obtained  from  the  hybrid 
algorithm. 

The  theory  has  been  developed  and  is  reported  in  Section  III.  Numerical 
implementation  has  been  initiated  and  will  be  continued  under  a  separate 
Task. 

3.  In  addition  to  performing  the  technical  tasks  described  under  1  and  2, 
substantial  effort  has  been  exerted  on  the  preparation  of  manuscripts 
suitable  for  publication.  The  material  in  Secs.  II  and  III  has  been  written 
up  in  such  a  manner  as  to  suit  this  objective.  Each  section  is  complete  in 
itself  and  any  references  within  the  text  refer  to  that  section  only. 
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Section  II 

Spectral  Options  to  Improve  the  Paraxial  Narrow  Gaussian  Beam  Algorithm 
for  Critical  Reflection  and  Head  Waves 
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Spectral  Options  To  Improve  the  Paraxial  Narrow  Gaussian 
Beam  Algorithm  for  Critical  Reflection  and  Head  Waves 

XJ.  Gao  and  LB.  Felsen 


I.  Introduction 

Wave  propagation  over  long  distances  in  complex  surroundings  that  vary  slowly 
over  the  scale  of  the  local  wavelength  is  of  importance  in  many  disciplines  dealing  with 
the  terrestrial  environment.  Direct  numerical  methods  are  unsuitable  because  the 
problem  is  computationally  large.  Therefore,  analytical  models  based  on  various  types 
of  propagators  have  been  explored  to  render  the  problem  tractable  -  an  exercise  that 
requires  a  trade-off  between  analytical  approximations  and  numerical  complexity. 
Ray  methods  are  favored  in  propagation  algorithms  because  they  are  easily 
implemented  numerically,  but  they  suffer  deficiencies  in  ray  field  transition  regions 
associated  with  focusing  effects  (caustics  and  foci),  shadow  boundaries,  diffraction, 
critical  reflection,  etc.  These  defects  can  be  formally  removed  on  replacing  the  ray 
fields  with  paraxial  (PA)  Gaussian  beam  fields  (GB),  because  such  beams  do  not  fail 
in  transitional  domains,  and  tracking  them  numerically  is  only  slightly  more 
cumbersome  than  ray  tracking.  Moreover,  by  an  extension  of  real  ray  theory,  the  PA 
Gaussian  beam  can  be  modeled  as  a  bundle  of  complex  rays  emanating  from  a 
complex  source  point  (CSP),  thereby  permitting  its  tracking  through  a  simple 
extension  of  the  procedures  for  real  rays.  These  aspects  have  been  well  documented 
in  the  literature  [1].  Since  the  final  field  is  constructed  by  beam  shooting,  one  avoids 
the  need  for  eigenray  search  which  would  connect  the  CSP  and  the  observer  along  a 
direct  complex  ray  path.  Instead,  the  field  at  the  observer  is  synthesized  by  paraxial 
■  expansion  about  each  on-axis  beam  field  (each  beam  axis  follows  a  real  ray),  thereby 
simplifying  the  algorithm  substantially. 

The  major  task  in  employing  any  scenario  for  propagation  modeling  is  to  assess  its 
accuracy.  This  has  to  be  done  by  comparing  data  from  the  proposed  algorithm  with 
data  from  an  independently  generated  "reliable’’  solution,  preferably  a  "benchmark" 
with  known  accuracy.  The  simplest  benchmarks  are  draped  around  idealized 
canonical  configurations  that  incorporate  relevant  wave  phenomena  within  a  formally 
exact,  and  numerically  computable,  framework.  If  discrepancies  arise,  careful  study  of 
these  canonical  solutions  may  reveal  the  reasons  behind  the  shortcomings  of  the 
approximate  model,  and  ways  to  repair  the  deficiencies. 

In  the  PA-GB  stacking  scheme,  one  of  the  critical  beam  parameters  is  the  beam 
waist.  Wide-waisted  beams  (in  terms  of  wavelength)  have  good  PA  behavior  but  are 
difficult  to  propagate  through  a  complex  environment,  whereas  the  opposite  is  true  for 
narrow-waisted  beams.  Even  if  defects  can  be  tuned  out  initially  by  numerical 
experiment,  such  tuning  generally  will  not  remain  intact  as  the  beams  progress 
through  various  refraction,  reflection,  diffraction,  etc.,  encounters.  For  many 
applications,  especially  in  seismology  and  bottom  interacting  ocean  acoustics,  critical 
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reflection  and  generation  of  head  waves  play  an  important  role.  Narrow  PA-GB 
methods  are  known  to  fail  here  [2].  A  detailed  canonical  prototype  study  has  been 
performed  previously  to  identify  the  cause  of  the  failure:  narrow  PA-GB  are 
spectrally  deficient  in  their  ability  to  cover  both  the  specularly  reflected  and  head  wave 
contributions,  whereas  wide  PA-GB  are  not  [3],  However,  wide  PA-GB,  as  noted 
above,  are  numerically  inefficient  and  cannot  be  well  adapted  to  non-canonical 
propagation  events. 

In  the  present  study,  we  have  returned  to  the  PA  narrow  beam  critical  angle 
problem  for  a  line-source-excited  critically  reflecting  interface  separating  two 
homogeneous  half  spaces  (either  for  horizontally  polarized  shear  (SH)  waves  in 
elastic,' or  for  compressional  (P)  waves  in  fluid,  environments),  searching  for  ways  in 
which  to  flesh  out  the  deficient  spectra  without  sacrificing  the  PA  assumption.  As 
may  be  anticipated,  this  endeavor  turns  out  to  be  accompanied  only  with  partial 
success.  Initially,  the  investigation  has  been  concerned  with  assessing  the  following 
options:  1.  inclusion  of  the  lateral  (head)  wave  in  the  total  reflection  domain,  outside 
the  critical  reflection  transition  region,  where  lateral  wave  asymptotics  is  invalid  [4];2. 
inclusion  of  lateral  shifts  in  the  reflected  beam  fields,  thereby  yielding  in  the  critical 
reflection  domain  two  reflected  contributions  (one  launched  near  the  critical,  and  one 
near  the  specular,  angle).  3.  inclusion  of  the  lateral  wave  and  one  (the  near-specular) 
laterally  shifted  reflected  ray.  By  systematic  PA  asymptotics  as  well  as  justifiable 
phenomenological  considerations,  each  of  these  options  accounts  for  the  critical 
reflection  and  lateral  wave  effects  individually.  For  the  laterally  shifted  beams,  the 
interface  plane  wave  reflection  coefficient  is  incorporated  into  the  phase,  and  the 
resulting  phase  is  made  stationary  to  determine  the  beam  axis  trajectory.  Numerical 
integration  of  the  exact  spectral  integral  has  been  taken  as  the  the  reference  solution, 
with  which  these  options  are  compared. 

In  these  studies,  the  PA  was  retained  in  its  conventional  form  wherein  the  phase  at 
an  off-axis  observation  point  vith  respect  to  a  particular  beam  in  the  stack  is  evaluated 
by  perturbation  (up  to  second  order  in  the  off-axis  distance)  about  the  on-axis  field. 
Tlie  results  still  showed  deficiencies,  which  were  felt  to  be  attributable,  at  least  in  part, 
to  the  conventional  PA  scheme,  which  is  strained  beyond  its  limits  for  narrow-waisted 
beams.  To  improve  the  approximation  without  sacrificing  beam  shooting,  the  off-axis 
results  were  recalculated  by  retaining  the  exact  complex  distance  from  the  C$p  to  the 
observer.  This  modification,  which  can  easily  be  incorporated  into  the  algor  \nm,  has 
been  found  to  reduce  the  discrepancies  with  the  reference  solution,  and  it  is  therefore 
recommended.  . 

v 

In  what  follows,  the  spectral  integral  reference  solution  is  stated  in  Sec.  II.  Section 
III  summarizes  the  PA-GB  algorithm,  both  in  its  conventional  (paraxial  CSP)  and 
improved  (full  CSP)  versions.  The  augmentations  of  PA-GB  described  by  options  1.- 
3.  above  are  treated  in  Sec.  IV,  and  the  corresponding  numerical  comparisons  for 
typical  samples  extracted  from  an  extensive  set  of  data  are  discussed  in  Sec.  V. 
Concluding  remarks  are  found  in  Sec.  VI. 
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II.  Spectral  Integral  Reference  Solution 


We  consider  acoustic  wave  propagation  in  the  presence  of  a  plane  interface  (z=0) 
that  separates  an  upper  half  space  (z  >  0;  medium  1)  from  a  lower  half  space  (z  <  0 ; 
medium  2)  (Fig.  1).  Time-harmonic  excitation  (with  suppressed  dependence  exp(-iwt) 
is  from  a  y-directed  line  source  located  at  (xs,  zs).  By  plane  wave  spectral 
decomposition,  at  an  observation  point  (x,z)  in  medium  1,  the  total  field  due  to  the 
interface  is  given  by  [5], 


I  = 


-1 

47ri 


/ 

-OO 


m. 

«i(0 


M.0  (ft 


(i) 


where  the  integration  path  is  shown  in  Fig.  2(a),  £  is  the  spectral  parameter 
(wavenumber  along  x)  in  the  plane  wave  continuum, 


•  V<£)  =  £(x-xs)  +  «i(z  +  zs),  Ki(0  =  \/k? -£2,  Im (/q) >0,  i  =  1,2  (la) 


and 

...  =  A1/c1(0-A2/c2(0 
Ai«i(C)  +  A2iC2(0 

is  the  interface  reflection  coefficient.  The  wave  propagation  speeds  in  medium  1  and 
2  are  v^d  v2,  respectively,  the  corresponding  wavenumbers  are  k1>2  =  w/v1)2,  and 
the  densities  are  p1>2.  An  angular  spectrum  reprsentation  in  terms  of  the  angle  9  (Fig. 
2(b))  is  obtained  via  the  mapping 

£  =  ki  sin0  (2) 


Ai  =  pjV2,  i  =  1,2  for  SH  waves  (elastic)  (lb) 
Ai  =  p2,A2  =  pi  for  P  waves  (fluid)  (ic) 


The  integration  contours  in  the  complex  £  and  9  planes,  consonant  with  the  radiation 
condition,  are  shown  in  Fig.  2. 


III.  The  Gaussian  Beam  Method 

The  conventional  GB  method  [1],  whose  formulations  have  been  critically 
examined  in  [5]  and  [6],  is  briefly  summarized  here.  By  that  method,  the  incident  field 
from  the  line  source  is  approximated  by  a  continuum  of  PA-GB  whose  amplitudes  are 
adjusted  so  as  to  reproduce  the  asymptotic  far  field  from  the  exact  angular  spectrum 
(see  (1),  with  (2),  and  F  replaced  by  1).  Tne  angular  spectrum  PA-GB  integral  is  i’nen 
discretized  for  numerical  evaluation,  with  Aa  denoting  the  angular  separation  between 
adjacent  beams,  which  are  assigned  a  beam  waist  w0  =  (2b/k1)1/2,  with  b  denoting  a 
parameter  equivalent  to  the  Fresnel  length,  b  and  Aa  are  at  the  disposal  of  the  user, 
and  this  is  at  the  root  of  the  arbitrariness  that  besets  this  method  [6].  Considering 
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propagation  to  the  right  as  in  Fig.  3(a),  the  angular  spectrum  interval  may  be 
truncated  at  and  cS2\  assuming  (subject  to  numerical  test)  that  beams  outside 
that  interval  make  a  negligible  contribution.  This  leads  to  the  following  PA-GB 
expansion  of  the  integral  (1): 

I-  E  GnA a  ,  -£<ai  < a2  <£,  < <  1  (3) 

a0)  *  ^ 

Each  incident  and  reflected-GB  can  be  modeled  by  a  source  at  a  complex  location 
(CSP),  with  the  complex  displacement  for  a  given  beam  waist  specified  by  the 
paramtti  \  Since  the  asymptotic  reflected  field  derived  from  the  spectral  integral  in 
(1)  is  equivalent  to  the  field  emanating  from  an  image  source  located  as  in  Fig.  3(b)  in 
the  infinitely  extended  medium  1,  with  the  strength  of  the  image  source  given  by  the 
spectral  reflection  coefficient  T,  one  may  write  the  M  CSP  asymptotic  approximation 
for  each  reflected  beam  as  follows  [5]: 

\Zki(b-iS) 

$n=$(on)  =  $(0)  = - 2it . -  exp(ik1S-k1b)  (3a) 

l 

(3b) 


where  Rn  is  the  complex  radial  distance  from  the  n-th  beam  waist  center  to  the 
observer  at  (x,z), 

Rn  =  V  (x-Xsn)2  +  (z  +  Z^)2 ,  Re(R„)  >  0  (3c) 

and  0n  denotes  the  complex  observation  angle  given  by 

x-xsn 

sin0n  =  -f-  (3d) 

Rn 

The  beams  are  launched  from  the  initial  surface  with  radius  S  in  the  "far  zone"  of  the 
line  source  located  at  (xs,zs).  The  n-th  beam  waist  center  on  s  is  located  at  (xsn,zsn), 
where  (see  Fig.  3(a)) 

xsn  =  xs  +  S  sin  an,  zsn  =  zs  +  Scosan  ^  (3e) 

with  an  denoting  angular  location.  In  tarns  of  these  quantities,  the  corresponding 
CSP  (xsn,zsn)  for  the  n-th  beam  is  given  by 

Xsn  =  xsn  +  ibsinan,  z^  =  zsn  +  ibcoso:n  (3f) 

The  complex  coordinates  are  schematized  in  Fig.  3. 


GnH^msinSn) 


2 IT 


ktRn 


exp 


ik.tRn 


nr 
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The  conventional  PA  is  recovered  from  the  full  CSP  result  in  (3b)  by  the  following 
approximations: 


1.  expanding  Rn  in  the  phase  to  second  order  in  the  off-axis  distance  dn: 

1  dl 


Rn~^an  +  9  _  -u 

2.  Ton  *  lb 


2.  expanding  Rn  in  the  amplitude  factor  to  first  order: 


~^ar, 


(4a) 


(4b) 


3.  replacing  the  complex  angle  0n  by  the  real  beam  angle  an  in  the  reflection 
coefficient  P: 


T(ki  sin^n)~  P(ki  skflO  (4c) 

With  reference  to  Fig.  3(a),  one  identifies  the  various  geometrical  quantities:  dn  is  the 
perpendicular  displacement  from  the  point  (x^.z^)  on  the  beam  axis  along  an  to  the 
observer  at  (x,z),  while  r  is  the  radial  distance  from  the  image  source  to  (x,z)  and  6  is 
its  angular  displacement  from  the  z-axis.  Explicit  expressions  for,  and  relations 
between,  these  quantities  are: 


r  =  V^lo  +  do  ,  =  «o  +  tan-1 

do_ 

ra0 
k  ) 

■  (4d) 

where  do  >0  if  x>xa0,  do  <0  if  x  <  xa0 

dn  =  rstn^on) 

(4e) 

ran  =  rCOs(5-^)-S 

(4f) 

xan  —  Xj  +  r  sin  otn  ,  zan  —  zs  4* 

(4g) 

Thus,  the  complex  rad.*u  distance  to  (xan?zan)  is 

^•an  y/ (*an  "  *sn )  ^  (^an  d"  2sn)^ ,  Re(Ran)  >  0  (4h) 


Implementation  of  the  algorithm  in.  (3)  via  the  conventional  PA-GB  model  in 
4(a)-4(h),  leads  to  the  difficulties  discussed  in  [5].  For  narrow-waisted  beams, 
substantial  discrepancies  exist  between  the  PA-GB  and  reference  solutions,  with  the 
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former  missing  essentially  the  (reflected  wave)  -  (head  wave)  interference.  This  can 
be  attributed  to  the  inability  of  the  narrow  PA-GB  spectrum  to  account  for  both  wave 
phenomena.  It  is  therefore  suggested  to  augment  the  algorithm  in  (3)  by  including 
each  of  these  phenomena  individually.  The  options  to  be  pursued  have  been  noted  in 
Section  I.  Subsequently,  we  shall  show  that  the  full  CSP-GB  reduces  the 
'discrepancies. 

IV.  Augmented  Conventional  PA-GB  Method 

A.  PA-GB  plus  head  wave 

In  this  option,  we  include  the  head  wave  contribution  excited  by  the  n-th 
conventional  PA-GB  in  (3),  with  (4).  The  domain  of  existence  of  the  n-th  head  wave 
can  be  ascertained  by  returning  to  the  exact  spectral  integral  for  the  line  source  in  (1), 
changing  the  line  source  input  into  the  n-th  GB  input  through  replacement  of  the  real 
source  coordinate  (xs,Zs)  by  the  n-th  beam  complex  source  (CSP)  coordinate  (xsn.Zsn)> 
and  performing  an  asymptotic  evaluation  of  the  resulting  exact  integral  for  the  total 
field  excited  by  the  n-th  incident  GB.  This  requires  deformation  of  the  original 
integration  path  into  the  steepest  descent  path  (SDP)  through  the  complex  saddle 
point  £n  (Fig.  2).  If  the  branch  point  at  £b  =  ^2  is  intercepted  during  this  deformation, 
(this  occurs,  roughly  speaking,  for  observation  angles  6>0C,  where  0C  =  sin"1  (vj/v 2)  is 
the  angle  of  critical  reflection),  a  branch  cut  contribution  Ibn  must  be  added  to  the 
SDP  contribution;  asymptotically,  the  former  generates  the  reflected  GB  Isn  whereas 
the  latter  generates  the  head  wave  excited  by  that  GB.  Asymptotics  are  valid  when 
the  branch  point  and  saddle  point  are  sufficiently  well  separated.  A  measure  of  the 
separation  is  provided  by  the  numerical  distance  ND,  and  it  has  been  found  that 
ND>2  is  adequate. 

A  detailed  discussion  is  given  in  [5].  Correcting  the  omission  of  the  density  and 
velocity  ratio  in  the  SH  wave  formula  of  that  paper  (the  correct  formula  is  given  in 
(lb)),  it  is  found  that  the  contribution  excited  by  the  n-th  incident  beam  is  given  by: 


In  “  I$n  +  Ibn  U  (£r  -  k2) 

(5) 

fO  X<0 

U  (x)  =  4  1/2  X  =  0 

(5a) 

U  x  >  0 

where  £r  denotes  the  intersection  of  the  SDP  with  the  real  axis  in  the  complex  £-plane, 
and  Isn  —  Gn  is  given  in  (3b).  The  branch  point  contribution  is  given-by 


1  A2  y/  tan0c  ik,R.cosAs-i-j- 
A!  (iqRn  sinAn)3/2 


(5b) 
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The  following  definitions  have  been  used  here, 


sinAn  =  cosflc  -  --  -^Zsn-sin^e 
Rn  Rn 


X-X, 


cos  An  =  — — sin0c  + 
Rn 


Z  +  Zsn 
Rn 


COS  dr 


Re[k1RnsinAn]3/2>0 

The  numerical  distance  is  given  by  [7], 

ND  =  \{\ 

where 


?  =  i(-2ikn) 2  A 
27 


/?  = 


sink 


z  +  z, 


sn 


Rn 


•COSk 


X-X, 


sn 


Rn 


/  sin20c 


7"T 


Z  +  Zsn 


X-X, 


+ 


sn 


Rncos 39c  Rnsin30c 


0C  =  sin'1 


21 

v? 

\  j 


and  Rn  is  given  by  (3c). 


(5c) 


(5d) 

(5e) 


(5f) 

(5g) 

(5h) 


(5j) 


B.  Laterally  shifted  PA-GB 

This  option  is  based  on  an  alternative  procedure  for  the  asymptotic  evaluation  of 
the  spectral  integral  (1).  Instead  of  treating  the  reflection  coefficient  T  as  a  spectral 
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amplitude  as  in  Sec.  IV  A,  it  is  now  incorporated  into  the  phase  via  the  identity 
r  =  exp((n  r),  [8].  This  alters  the  stationary  phase  (saddle  point)  condition  and  yields 
modified  reflected  ray  trajectories  that  contain  an  additional  path  segment. 

For  a  ray  incident  beyond  the  critical  angle  ( 0>9C ),  /nT  is  imaginary  and  the 
additional  segment  can  be  interpreted  as  a  real  lateral  shift  along  the  interface  (see 
Fig.  4).  For  incidence  below  critical  (9<9C),  /nT  is  real  and  the  lateral  shift  is 
complex;  in  this  domain,  it  is  physically  more  transparent  to  retain  T  as  an  amplitude 
as  in  Sec.  IV  A.  It  turns  out  that,  two  laterally  shifted  paths  exist  beyond  critical 
incidence. 

Pursuing  the  strategy  outlined  above  for  the  n-th  GB  in  the  beam  stacking  scheme, 
we  write  (1)  in  the  CSP  form 

1  =  "XT  /  — brexP(iki'if(£)}<£  (Q 

*i(0 


where 


m«) = m  -  ¥ r  (o 

(6a) 

<K0  -  f(X-Xs)  +  Kl(z  +  Zs) 

(6b) 

A(()  =/nr«) 

(6c) 

Introducing  a  normalized  spectral  parameter 


(7a) 


and  determining  the  saddle  points  rnj  for  6  >  9C  from 


—  ¥(rs)  =  (x-xs)- 


v^T 


(z  +  zs)  -  A  =  0 


where 


(7b) 


_i_  _1_  dT 
ki  T  dr 


T  *  T, 


(7c) 


is  the  lateral  shift,  one  obtains  from  the  conventional  saddle  point  technique  [9]  the 
asymptotic  result 
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here, 


I|n«- 


-2  7T 


4 

2 


1  _  e«ki«(rN) 


V^nj 


^Onj)  =  ['(X-Xnj  +  VT^Cz  +  ZnsH/nrCr)]^ 


(8) 


(8a) 


d2tf(r) 

_  1  (-  ,  z  \ 

i 

S 

1  uT 

2  i  d2r‘ 

dr2 

(1.r2)3/2  +  ?ns) 

r«r 

"ki 

r«i 

r  dr 

r  d*2 

(8b) 


In  keeping  with  the  beam  shooting  strategy  that  avoids  solving  the  full  saddle  point 
condition  in  (7b)  for  a  given  observation  point  (x,z),  we  invoke  the  paraxial 
approximation.  If 

£<0)  =  ki  sin  On  (9a) 


is  the  saddle  point  for  an  observer  located  at  (xa,Za)  on  the  axis  of  the  laterally  shifted 
beam  which  is  incident  at  the  angle  an,  then  the  saddle  point  for  an.  off-axis  observer 
at  (x,z)  is 

(  =  ( f>  +  l  (9b) 

The  on-axis  and  off-axis  observation  points  are  related  as  follows 

x  =  Xar  +  dn  cosan  ,  z  •=  zan  -  dn  sina„  (9d) 


and 


dn  =  rsin(0  -  £*„)  -A  cosan 


where  r  and  0  are  given  by  (4d).  The  paraxial  expansion  of  the  phase  function  $(£;x,z) 
in  (6a)  into  power  series  about  (fj^x.z)  involves  up  to  quadratic  terms  in  the 
parameters  dn  and  S, 


¥(&X,  z)=^(^(n0);xan,zan)  +  ^~^(d0);x: 


'an> 


n. 

-^an)  4"  “’^(^bxan,  Zan)6 


2  d"  <9dn 


+  ^)2  ^Cri0):Xan4an  ) 


(10) 
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It  is  easily  verified  that 


3  v  ^  & 

3dn 


'®r(£n^>*an»zan)  —  _,2  ^(£a>x a»za)  — 0 


ad 


36 


*  (fSVan.Zan)  =  0  (saddle  point  condition) 


■  -dr- 

3an36  C.QSCCn 


(10a) 


(10b) 


(10c) 


kl5*  ^  (^n^JXan.Zan)  =  - 
do 


1 _ (z  +i  N 

3  Vzan  +  zsn;  ‘  2 


d2 


kicos^a 


=  Jli£|  +ld2r 


d£z 


(lOd) 


d0) 


de2  ^r(e)  "  [r  dej  ’  r  de2 


(10e) 


The  saddle  point  equation  in  (10b)  generally  has  two  solutions,  but  only  one  of  these 
turns  out  to  be  important  for  each  beam  element.  Inserting  the  above  expressions  into 
(10)  yields 

ki«'(£:x,z)  *  kl¥(£h0);Xan:zan)  +  dn5 

cosan 


1_ 

2 


zan  zsn 
kiCOS3^ 


+  1 


d20r(€<°>) 


<r 


(ii) 


Since  ^  +6  is  the  saddle  point  of  the  function  ^(f  ;x,z ),  it  must  satisfy  the  saddle 
point  condition 


kla»  *(£;x’Z)=:  o, 


which  furnishes  the  solution  for  6, 
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6  = 


cosa 


Z*n  +  Zsn  .  .  &&(£$) 


1*1 


kiCOS3an 


+i 


ae 


(12) 


Thus,  from  (11)  and  (12), 


k1'£(£;x,z)  =  k1*(^S9;xa,za)  +  j 


Zan  +  Zsn  .  . 

- 5 -  +1 

kiCOSOn 


de2 


(13) 


or  altematjV  \y  (see  (7a)  and  (13)) 

k!^(rpc,z)  =  k^^pc.z)  + 


(Mn)2 


kiRan  +  icos2an-d 


di2 


(14) 


where  Ran  is  given  in  (4h),  and 

d2V>r(rn) 

dr2 


■ 

f  > 

i  dr 

2  .  i  d2r 

r  dr 
\  > 

r  dr 2 

Also 


j£ 

dr2 


ki-TT^nPC.z)  =  * 


kiRan  +  icOSOtn 


d20r(^n) 


dr2 


/ 


C0S2On 


(15) 


Insertion  of  (14)  and  (15)  into  (8)  completes  the  paraxial  treatment  of  the  laterally 
shifted  beams  for  0>dc.  For  6<6C  one  retains  the  PA-GB  in  (3b)  with  the 
approximations  in  (4a-c).  The  beam  sum  in  (3)  comprises  all  contributions 
(nonshifted  for  &<9C  and  shifted  for  9>BC)  that  fit  into  the  interval  <an  <oP\ 

C.  Reflected  laterally  shifted  PA-GB  plus  head  wave 

This  option  is  a  hybrid  combination  of  those  in  Secs.  IVA  and  IVB,  motivated  by 
the  observation  that  the  laterally  shifted  ray  incident  near  the  critical  angle  does  not 
account  for  the  head  wave  as  well  as  the  branch  point  contribution  in  the  non-shifted 
formulation,  whereas  the  laterally  shifted  ray  incident  near  the  specular  reflection 
angle  has  little  effect  on  the  head  wave  but  does  improve  the  totally  reflected 
contribution  (the  reflection  coefficient  T  is  not  necessarily  slowly  varying  (i.e.  an 
amplitude  function)  for  strongly  oblique  incidence).  As  before,  incident  ray  fields 
ahead  of  and  near  the  critical  angle  are  modeled  by  PA-GB. 
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The  GB  implementation  of  the  integral  (1)  here  was  initially  structured  as  follows. 
The  observational  domain  O<0<y  was  divided  into  three  subintervals:  l:O<0<#i; 

2:Qi<Q<&i;  3:^<0<y.  In  interval  1,  the  nonshifted  PA-GB  alone  were  used,  in 

interval  2,  the  PA-GB  were  augmented  by  the  (asymptotically  approximated)  head 
wave,  while  interval  3  was  covered  by  the  laterally  shifted  PA-GB  as  in  (8)  plus  the 
(asymptotically  approximated)  head  wave.  Accordingly, 


In=<! 


sn  > 


9<6<6l 


Isn  +  Ibn  >  01<0<02 


(16a) 

(16b) 


I  In  + 


bn 


'h  <6< 


7 r 


(16c) 


where  Isn  and  Ibn  are  given  by  (3b)  and  (5b),  respectively,  while  is  given  by  (8). 
The  angle  6l  >  6C  corresponds  to  a  numerical  |  ND  |  =2,  whereas  the  angle  £>  >  0C 
corresponds  to  a  numerical  distance  |  ND  |  >2  which  is  problem  dependent.  From 
numerical  experience  for  various  ranges  of  medium  parameters  and  source  locations, 
we  chose  1.02 9C. 

The  choice  of  the  three  intervals  above  was  regarded  to  be  in  accord  with  the 
phenomenological  spirit  of  the  hybrid  formulation.  Nonshifted  PA-GB  plus  branch 
point  asymptotics  can  account  for  phenomena  near  the  critical  angle  because  the 
branch  point  (head  wave)  and  saddle  point  (reflected  beam)  in  the  spectral  integrand 
are  so  close  to  one  another  that  the  spectral  coverage  of  both  overlaps.  As  9  increases 
away  from  9C,  the  two  contributions  become  distinct  and  are  known  to  be  poorly 
accomodated  by  narrow  nonshifted  PA-GB.  This  deficiency  is  to  be  overcome  by  the 
laterally  shifted  reflected  PA-GB  for  9>&2,  where  $2  is  only  loosely  defined.  When 
this  scenario  was  implemented  numerically,  tests  were  run  as  well  wherein  the  interval 
2  in  (16b)  was  eliminated,  and  interval  3  in  (16c)  extended  from  9v<9<-nj2.  The 
results  obtained  so  far  were  practically  unchanged,  thereby  indicating  that  for  the 
problem  parameters  explored  here,  the  laterally  shifted  PA-GB,  which  accounts 
essentially  for  specular  reflection  away  from  the  critical  angle  can  be  safely  extended 
up  to  the  |  ND  |  =  2  limit.  This  simplifies  the  algorithm,  which  thus  comprises  (16a) 
and  (16c),  with  corresponding  to  |  ND  j  =  2. 


V.  Numerical  Results 

Extensive  numerical  results  have  been  generated  for  the  analytical  models  in  Secs. 
Ill  and  IV.  when  applied  to  typical  elastic  rock  interface  or  to  ocean-(fluid-bottom) 
environments.  Selected  sets  of  data  are  shown  in  Fig.  5,  with  primary  emphasis  on 
how  the  results  are  affected  by  the  choice  of  beam  widths  in  the  Gaussian  stack.  The 
half-width  w0  of  the  beam  waist  at  a  level  where  the  amplitude  equals  (1/e)  = 
1/2.71828  times  the  maximum  amplitude,  is  related  to  the  beam  parameter  b,  which 
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measures  the  Fresnel  distance  corresponding  to  wo,  as  follows: 

b  =  7TW§/A! 


(17) 


where  A^  is  the  wavelength  in  the  source  medium  1.  The  angular  spacing  Aa  between 
beams,  for  specified  w0  or  b,  is  such  that  adjacent  beams  overlap  at  their  (1/e) 
amplitude  points.  The  total  number  of  beams  in  the  stack  is  determined  by  any  one  of 
the  following  three  criteria:  a)  |  6 -ecu  I  >/ tt/2,  where  an  is  the  beam  axis  angle  and  6 
is  the  angle  of  the  ray  that  connects  source  and  observer  (the  n=0  beam  has  angle 
ao  =  0);  b)  (amplitude  of  beam  with  on  + 1)  <  (10"4  times  n  =  0  to  N  beam  sum);  c) 
on  >tt/2.  Note  that  N  may  be  different  for  aR  <  0(N  =  Ni)  and  an  >  0(N  =  N2),  so 
that  the  total  number  of  beams  is  Ni+N2+1).  Each  beam  is  treated  by  the  CSP 
method,  which  involves  the  complex  image  point  (xsn,  -zsn),  with  xsn  =  xsn+  bsinan, 
zsn  =  zsn  +  ibcosa„  (this  is  a  convenience,  not  a  restriction;  alternatively,  the  beams 
could  have  been  tracked  from  the  source  to  the  observer  via  the  interface);  here, 
(xsn»zsn)  locates  the  n-th  beam  waist  center.  Other  relevant  quantities  have  been 
defined  in  (3e). 

We  recall  that  in  the  PA-GB  option,  the  beams  alone  (either  via  the  conventional 
or  the  improved  PA)  are  included  (see  (3)  and  (4)).  In  the  PA-GB  plus  head  wave 
option,  the  PA-GB  are  included  as  above  (see  (3)  and  (4)),  and  the  head  wave  in  (5b) 
is  added  to  the  n-th  beam  when  each  of  the  following  two  criteria  is  satisfied:  a)  9  >  9C, 
where  6G  is  the  critical  angle;  b)  (ND)n>2.0,  where  (ND)„  is  the  numerical  distance  in 
(5f).  While  (ND)n>2  has  been  found  safe  for  validity  of  branch  point  asymptotics  over . 
a  broad  parameter  range,  one  may  "tune  up"  the  restriction  to  lower  ND  values  by 
numerical  experiment  in  specific  cases. 

When  the  beam  parameter  b  is  larger  than  50  Aj,  the  numerical  results  become 
unstable  and  indicate  that  the  angular  interval  Aa  must  be  reduced  below  the  (-1/e) 
overlap  criterion  by  the  addition  of  more  beams.  The  difficulty  is  caused  by  the 
increasingly  rapid  variation  of  the  phase  term  UqRn  in  (3b)  with  an,  as  b  increases. 
For  example,  when  500Ai  <b  <  1000A1}  it  has  been  found  that  Aa= 0.004  is  needed  for 
stability  whereas  the  1/e  overlap  criteria  gives  Aa~  0.0178.  It  may  also  be  noted  that 
the  laterally  shifted  beams  near  the  critical  angle  require  tighter  sampling  than  the 
nonshifted  PA-GB.  Again,  when  500Ax  <  b  <  lOOOAx,  stability  sets  in  when 
Aa=0.00178.  As  before,  the  difficulty  arises  from  the  very  rapidly  fluctuating  phase 
when  the  reflection  coefficient,  which  varies  rapidly  near  the  critical  angle,  is 
incorporated  into  the  phase.  The  sensitive  regions  are  near  the  critical  angle  and  near 
the  zero  reflection  (Brewster)  angle  (when  that  exists);  when  sampling,  is  the 
synthesized  results  are  not  smooth  in  these  regions. 


Addressing  the  above  difficulties  by  increasing  the  number  of  beams  led  to  the 

For  the  conventional  nonshifted  PA-Gp- 
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interval  was  determined  as  follows: 
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Ac  ~  0.2 


7rb 


1 1/2 


Aa  =  0.15 


Ai_ 

flrb 


1/2 


50Ai  <  b  <  50QAi 


500A!  <b<500QA« 


The  corresponding  number  of  beams  is 

b  =  10Q\i(w0  =  5.6AO,  Ni  +  N2  ~  50 
-  b  =  IOOOA^wq  =  17.8A!),  Nx  +  N2  « 130 
b  =  5000A!(wo  =  39.9AO,  Ni  +  N2  ~  580 


(17a) 


(17b) 


(18) 


For  the  nonshifted  PA-GB  plus  head  wave  option,  Aa  and  Ni  +  N2  are  the  same  as 
above.  For  the  hybrid  option,  more  beams  are  required  for  numerical  stability.  By 
numerical  tests,  we  chose 


Aa  =  0.06 


7rb 


b  <  1000A! 


(19a) 


Aa  “  0.04 


1000 Ax  <b<5000A! 


with  the  corresponding  number  of  beams 

b  =  100A1,  Nx  +  N2  =280 
•  b  =  IOOOAl  Ni  +  N2~890 
b  =  5000Ab  +  N2~980 


(19b) 


(20) 


Truncation  of  the  beam  staci  s  according  to  the  criteria  listed  after  (17)  must  be 
performed  with  care  when  the  bottom  medium  allows  for  Brewster  angle  effects 
where  the  plane  wave  reflection  coefficient  vanishes.  Although  the  reflected  beams 
very  near  the  Brewster  angle  0  =  &b  make  a  negligible  contribution  at  the  observation 
point,  the  beam  stack  may  be  truncated  there  only  when  these  contributions  remain 
negligible  in  a  large  enough  neighborhood  of  4-  We  have  found  that  such  a 
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neighborhood  is  adequately  defined  by  ^  ±  A,  where  A =0.01^. 

For  the  laterally  shifted  PA-GB  option  in  Section  IVB,  similar  caution  must  be 
exercised  concerning  truncation.  The  laterally  shifted  PA-GB  stack  yields  strongly 
contributing  clusters  near  the  critical  angle  6C  and  near  the  specular  reflection  angle  0 . 
Truncation  of  the  total  beam  stack  by  applying  the  above-cited  criteria  can  lead  to  the 
exclusion  of  the  other  cluster.  Therefore,  no  truncation  should  be  undertaken  in  the 
angular  interval  between  9C  and  6. 
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Fig.  1  Physical  configuration  of  line-source-excited 
half  space  with  interface  at  z=C 
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1  Conventional  and  Augmented  PA-GB 


Fig.  5,  SH  field  synthesis  by  conventional  PA-GB  and  its  augumented  options  for 
lint  source  in  the  presence  of  a  plane  interface  separating  two  elastic  half 
spaces 


Parameters: 

frequency  f*lhz 

wave  length  in  medium  1  A(  *  1 

density^  *  density/),  *  1 

Vi  ”  £V;/«c 

vj  *  11S2jv[ 

critical  angle:  te  *  27.6  Degs.  (x  &  1043  At) 
observation  plane:  ?»100  A. 

tt.........  ~  A  -  W.AAV 

numerical  distance  NT)  »2 


—  •  Reference  solution 
A  :  PA-GB  only 
Q  :  PA-GB  plus  head  wave 
+  :  PA-G3  with  lateral  shift 


(a) b“100A[  (w0»5.6Aj) 

(b) b-200Ai  (w0«7.98A,) 

(c)  b  -  500  At  (w0  -  12.6A.) 

(d) b*1000A,  (wo-17.8A,) 

(e)  b  *2000  A,  (w0  »  2S.2A  t) 

(f)  b  *5000  At  (wj*39.9A|) 
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Conventional  and  Augmented  PA-GB 


Fig*  5  e) 
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Section  III 

Hybrid  Ray-Mode  Formulation  of  High  Frequency  Propagation 
in  a  Bilinear  Tropospheric  Surface  Duct 
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Hybrid  Ray- Mode  Formulation  of  High  Frequency  Propagation 
in  a  Bilinear  Tropospheric  Surface  Duct 

T.  Ishihara  and  L.B.  Felsen 

I.  Introduction 

Analytical  modeling  of  wave  propagation  In  the  earth’s  environment  poses 
formidable  problems  because  of  the  complexity  of  the  propagation  channel. 
The  troposphere,  which  provides  such  a  channel  for  electromagnetic  waves,  can 
generally  be  described  by  a  permittivity  e(>)  that  varies  vertically  as  well  as 
laterally  with  respect  to  the  earth’s  surface.  Of  special  Interest  are  conditions 
where  the  vertical  profile  of  e(r)  Is  such  as  to  permit  wave  ducting  adjacent  to 
the  earth’s  surface  or  at  moderate  heights,  with  lateral  variations  occurring  on  a 
scale  that  Is  gradual  compared  with  those  In  height.  Because  the  problem  scale 
Is  large  at  high  frequencies  and  for  long  propagation  distances,  direct  numerical 
modeling' Is  either  Inefficient  or  not  feasible.  Therefore,  an  effective  algorithm 
must  be  parametrized  In  terms  of  wave  objects  that  can  negotiate  long  range 
propagation.  Possible  choices  Include  ray  fields,  mode  fields,  parabolic 
propagators,  beams,  etc.,  either  for  single  species  models  or,  more  effectively, 
in  hybrid  combinations  that  seek  to  take  advantage  of  the  most  favorable 
features  of  each.  The  analytical  models  must  be  self-consistent  approximations 
of  the  Intractable  rigorously  formulated  problem.  To  ascertain  their  accuracy, 
they  must  be  tested  on  tractable  simpler  canonical  environments  that  allow 
generation  of  numerical  reference  solutions.  Moreover,  the  test  environments 
should  be  such  as  to  deviate  only  weakly  from  the  realistic  environment, 
thereby  validating  analytic  techniques  of  adiabatic  .adaptability  of  each  wave 
object  from  the  test  problem  to  the  realistic  case  without,  or  with  only  weak, 
coupling  to  other  wave  objects. 

The  present  study  seeks  to  apply  these  concepts  to  high  frequency 
electromagnetic  wave  propagation  in  a  tropospheric  surface  or  elevated  duct, 
Initially  without,  and  eventually  with,  lateral  inhomogeneities.  In  the  test 
environment,  the  permittivity  Is  assumed  to  depend  on  height  only,  and  the 
earth’s  Influence  Is  modeled  by  a  constant  surface"  Impedance,  thereby 
rendering  the  propagation  problem  separable  In  a  spherical  coordinate  system. 
Because  the  (vertical  eiectrlc  aipoie)  source  ana  observer  locations  of  lute  ret 
are  at  levels  above  the  earth  which  are  small  compared  to  the  earth's  radius,  it 
is  possible  to  Invoke  approximations  that  map  the  spherical  geometry  Into  an 
equivalent  rectangular  geometry.  A  normalized  field  solution,  usually  referred 
to  as  the  attenuation  function,  |s  then  obtained  In  a  spectral  Integral  form  that 
serves  as  the  starting  point  for  developing  desirable  parametrizations.  The 


35 


Science  and  Technology  Corp. 
Technical  Report  3021 


considerations  pertaining  to  this  strategy  have  been  well  documented  in  the 
technical  literature  [],  and  they  are  summarized  in  Sections  II  and  III.  In 
Section  IV,  detailed  consideration  is  given  to  a  bilinear  permittivity  profile 
which  has  also  received  much  attention  from  other  investigators  (].  After 
deriving  the  spectral  Integral  solution  for  that  special  case  (Section  rv.A), 
alternative  solutions  are  developed  in  the  form  of  an  expansion  over  normal 
(guided)  modes  propagating  parallel  to  the  earth’s  surface  (Section  IV.B),  a  ray 
field  expansion  (Section  IV.C),  and  a  hybrid  ray-mode  expansion  that  combines 
ray  fields  and  mode  fields  self-consistently  (Section  IV.D). 


II.  Formulation 


Adopting  a  spherical  r=(r,d,<f>)  coordinate  system,  we  consider  wave 
propagation  due  to  a  vertical  electric  dipole  located  near  the  earth’s  surface  in  a 
radially  varying  troposphere  modeled  by  a  permittivity  e(r).  The  electrical 
properties  on  the  earth’s  surface  at  r=  a  are  assumed  to  be  specified  by  a 
surface  impedance  Zs.  In  view  of  the  environmental  dependence  on  r  only,  the 
radial  dipole  source  excites  fields  that  are  azlmuthally  independent  about  Its 
axis.  If  the  polar  axis  of  the  spheric.al  coordinate  system  is  chosen  to  pass 
through  the  source  at  r=r/,  the  magnetic  and  electric  fields  are  specified  in  that 
system  by  the  components  and  Er, Eg,  respectively  (Fig.  1).  All  field 
components  can  be  derived  from  a  scalar  potential  U(r)=U(r,0)  via  the 
relations  [1,2] 


Er  = 


Eo  = 


- 1 

d 

,ln „  3U 

(i) 

rsin0 

dd 

sl"e  as 

1 

d_ 

, ,  au  i 

l ( J  in 

(2) 

e(r)r 

dr 

iw£(r) 

du 

dd 

(3) 

with  U(r)  determined  from  the  scalar  wave  equation 


d  i 
dr  e(r) 


3_ 

dr 


[c(r)(rU)J  + 


l 
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-rrSintf-rr-U 
d6  dd 
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k2(r)  =  w2e(r)/i0  -  w2e0/i0€s(r)  =  k2e3(r)  (4a) 

In  (4a),  e3(r)  Is  the  permittivity  normalized  with  respect  to  the  permittivity  e0 
in  vacuum,  and  k0  is  the  wavenumber  in  vacuum.  All  field  quantities  have  a 
suppressed  time  dependence  exp(-iwt).  The  Impedance  boundary  condition 

E*  =  -ZsH*  (5) 

implies  that 

-..L-  j-(f(r)rUl  =  lcue(r) ZSU  on  r=a  (6) 

e(r)r  Or 

At  |r  |->oo,  the  potential  must  satisfy  the  radiation  condition. 


Because  the  radial  locations  of  the  source  and  observation  points  satisfy  the 
Inequalities  (r//a)«l,  (r/a)<<l,  the  exact  equation  In  (4)  can  be 
approximated.  By  a  sequence  of  scalings  and  redefinitions  of  coordinates,  one 
may  transform  the  field  dependence  from  the  spherical  (r ,9)  coordinates  into 
an  equivalent  dependence  in  rectangular  (y,x)  coordinates  [1,2].  The 
approximate  source-excited  transformed  wave  equation  becomes 


a2  d 

dy2  +i  dx 


+  p(y) 


■0(x,y)  =  -<5(x)<5(y-y) 


(7) 


subject  to  the  boundary  condition 

8 


dy 


+  qip  =  0  at  y=0 


(7a) 


and  a  radiation  condition  at  infinity.  The  following  definitions  have  been 
utilized: 


x  =  m 6,  y  -  — (r-  a),  y*  =  —  (r'-a)  , 
m  m 


(8a) 


k  =  k(0), 


m  = 


q=lm[es(0)]1/2Zs/Zo 


(8b) 


The  "attenuation  function"^  is  derived  from  the  normalized  potential 

U/U0  =  2(7 rx)1/2  e^/-4  ,  (9) 


where  U0  is  the  potential  for  the  same  dipole  located  In  the  presence  of  a  plane 
perfectly  conducting  earth.  Because  this  normalization  extracts  the  "strong" 
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phase  variation  along  8->(x/m),  one  obtains  the  simplified  parabolic  equation  in 
(7).  The  equivalent  permittivity  function  p(y)  Is  given  by 


p(y)  =  y+7(y). 


l(y)  =  m2 


5(y) 

*,(o) 


1 


es(y)  =  cs(a+my/k)  (10) 


This  completes  the  formulation  of  the  problem. 


III.  Spectral  Integral  Solution 

The  reduced  eqution  In  (7)  for  the  attenuation  function  Is  separable  in  the 
(x,y)  coordinates.  The  most  general  form  of  the  solution,  from  which  various 
alternative  representations  can  be  derived,  Is  built  around  a  complex 
wavenumber  spectrum.  Introducing  a  complex  spectral  variable  (separation 
parameter)  t  (this  parameter  should  not  be  confused  with  the  deleted  time 
dependence  exp(-Icut)),  and  recognizing  that  the  (d/dx)  operator  Is  algebraized 
by  an  exponential  function  exp(lxt)  so  that  (l3/5x)— ►-t,  one  obtains  for  the  y- 
dependent  spectrum  the  one-dimensional  Green's  function  problem, 


!”£■  +  [p(y)-t]j  sy(yy;t)  =  -Ky-y*)  • 

.  (n) 

which  has  the  solution  [] 

,  .  ,  (r2(y<^)  +  R0fi(y<.«')]fl(y>,t) 

(12) 

where  the  Wronskian  W  is  given  by 

W(f„f2]  =  f^'-  f'f2  ,  r'  =  df/dt 

(12a) 

and  the  reflection  coefficient  R0  by 

f'(0,t)  +  qf2(0,t) 

R°  ”  f'(0,t)  +  qf^O,!) 

(12b) 

In  (12),  fi  and  f2  are  linearly  independent  solutions  of  the  so-urce-free  equation 
in  (11),  which  are  employed  in  the  combination  shown  to  satisfy  the  boundary 
conditions  aty-Q  and  y->oo,  respectively,  while  y<  and.  y>  denote  the  smaller 

or  larger  values  of  y  and  y',  respectively.  Spectral  synthesis 
following  solution, 

then  yields  the 
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=  ■  j— ,  §  e“*gy(y,y,;t)(lt  (13) 

—  2tt  i 

This  result  furnishes  the  starting  point  for  the  detailed  studies  that  follow. 

IV.  Bilinear  Permittivity  Profile 
A.  Spectral  integral 

We  shall  consider  for  the  equivalent  permittivity  p(y)  in  (10)  the  bilinear 
profile  (Fig.  2) 

7(y)  =  (l+p3)(yry)  for  0<y<y,  (14a) 

'y(y)  =  0  for  yf<y  (14b) 


where  n  is  a  parameter  to  be  determined  from  the  profile  expressed  in  the 
original  (r,0)  coordinates.  In  (14),  for  subsequent  convenience,  the  reference 
permittivity  has  been  placed  at  y=yjj  accordingly,  the  reference  permittivity 
<TS(0)  in  (10)  should  be  replaced1  by  e^yj).  If  we  write 


g3(y) 


1-bo"k"(y~y‘) 


y<yi 


(15) 


then  the  profile  parameters  In  (15)  and  (14)  are  related  via 


(16) 


To  reduce  the  equation  in  (11),  with  (14),  to  a  standard  form,  we  Introduce  the 
scaled  coordinate  and  spectral  variables  £  and  r, 


£  =  r+py,  r  -  (t-  (i+/z3)yj]//i2 


(17) 


from  which 

=  tr-y,  (17a) 

where  £=£,  corresponds  to  y=yj.  Then  the  generic  equation  for  the  source-free 
solutions  f=f12  in  (12)  becomes  the  Airy  equation 
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d£! 


-r-  er  =  0  T  <  Z  <  e, 


(18) 


Accordingly,  the  spectral  Integral  In  (13)  may  be  written  explicitly  In  terms 
of  the  Airy  functions  Al,  BI, 


■0(x,y) 


00  *f  (y<.t)'f(y>,t)  fxt 


—  J  — 

2/^1  -oo  l 


■R0(t)Ri(t) 


e,xMt 


(10) 


where 

f  (y<»t)  =  Al(r+/*y<)  +  R0Bi(r+^y<)  (19a) 

f  (y>»t)-=  Bl(r+^y>)  +  RiAi(r+/iy>)  (19b) 


with  the  spectral  reflection  coefficients  R0  and  Rj  that  arise  from  the  earth’s 
surface  at  y=  0  and  the  profile  slope  discontinuity  at  y=yif  respectively, 

_  /iAi/(r)+qAl(r)  ,  . 

^Bi'(r)+,BI(r)  U  ' 


wftt-y.) 

HBi'(T+ny{)  +  •■■■■■  ■  Bi(r+/iy,) 

Wjft-yi) 


w,'(t-y() 

P Al'( r  +/x yt) + Ai(r+/zy|) 

wntr'yiJ 


(19d) 


Alternative  representations  with  different  convergence  properties  will  now 
be  derived  from  the  spectral  continuum  In  (19). 

B.  Normal  Mode  Expansion 


The  Integrand  In  (19)  has  pole  singularities  tjn  at  the  zer's  of  the 
denominator  (the  Index  m  here  Is  not  to  be  confused  with  the  parameter  m  In 
(8b)), 

l-R0(tm)Ri(trn)  =  0  ,  m  =  1,2,3,---  (20) 


or  explicitly 

/zAl'(rm)+qAl(rm)  yl)Al,(rm+^yi)+w{(tm-yi)Ai(rm+;<y,) 

/iBl'( rm)  +qBI( r m)  “  fi wi(  tm-  yi)Bl/(rm+/iyj)+w((  t^-  y,)Bl(  r  m+^y,) 


By  deforming  the  integration  path  in  (10)  around  these  poles  In  the  upper  half- 
of  the  complex  t>-plane  (the  Integrand  converges  at  infinity  so  as  to  make  this 
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possible  [2])  and  invoking  the  residue  theorem,  one  generates  the  residue 
series 


^(x,y)  =  *7  E  - o - 

*  m=1  Ro(*ta)-^[Ro(Wt)]wi 

where 


(22) 


fm(y)  =  Ai(rm+/zy)  +  R0(tm)Bi(rm+^y)  (22a) 

Each  term  in  (22)  represents  a  mode  field  fm(y)  in  the  vertical  cross  section, 
which  propagates  along  x  with  propagation  coefficient  t^,.  The  profile  in  (14) 
creates  a  surface  duct  between  y=  0  and  y=y,  because  it  is  downward  refracting 
in  that  height  interval.  For  y>yit  upward  refraction  takes  place,  with 
consequent  leakage.  The  mode  set  tm  is  therefore  grouped  Into  surface  ducted 
(trapped)  modes  with  small  leakage  (Retm>p(yi)),  transitional  modes 
(Retm!5=Jp(yi)),  and  leaky  modes  (Retm<p(y())  (see  Fig.  3).  The  estimates 
relating  Retm  to  p(y,)  are  qualitative.  Losses  duj  to  leakage  out  of  the  duct  or 
due  to  a  dissipative  surface  Impedance  Z5  are  expressed  by  Im  t^,  which  is 
positive;’ this  controls  the  decay  rate  along  x  via  the  propagator  exp(ixtm)  in 
(22),  with  larger  values  of  Im  implying  stronger  damping. 

C.  Ray  Expansion 


The  normal  mode  expansion  in  Section  B  Is  not  always  the  most  convenient 
for  characterizing  high  frequency  wave  phenomena.  An  alternative  is  provided 
by  a 'ray  field  expansion  wherein  local  plane  waves  are  tracked  individually 
along  direct  and  multiple  reflected  trajectories  from  source  to  observer  [3,4]. 
Although  the  ray  expansion  can  be  phrased  rigorously  in  terms  of  generalized 
ray  spectral  Integrals,  the  principal  utility  of  a  ray  field  formulation  rests  In  the 
asymptotic  reduction  of  the  generalized  ray  integrals,  which  yields  the 
approximate  fields  that  satisfy  the  rules  of  ray  optics  in  Inhomogeneous  media. 
To  generate  the  generalized  ray  series  from  the  general  result  in  (10),  it  is  first 
necessary  to  decompose  the  oscillatory  Airy  functions  Ai  and  BI,  which  are 
natural  descriptors  of  the  (oscillatory)  normal  modes,  Into  traveling  wave  form. 
Next,  the  resonant  denominator  (1-RqR,)"1  in  the  integrand  of  (19)  is 
removed  by  power  series  expansion,  and  the  resulting  Integrand  is  rearranged 
into  a  sequence  of  combinations  of  traveling  wave  terms  that  form  the 
generalized  rays.  Finally,  asymptotic  reduction  produces  the  conventional  ray 
fields. 
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Implementing  the  above  scenario,  we  Introduce  the  traveling  wave  Airy 
functions 


Wj  2  =  Ai  IBI 

(23) 

and  thereby  rewrite  (19)  in  the  form 

• 

*t'(x,y<fy>)  =  -1  /  E(y<Jt)5iy>'t) 

4A*  c  l-Ro(t)R,(t) 

(24) 

where 

h  (y<,t)  =  w1(r+f/y<)  +  R0(t)w2(r+/iy<) 

(24a) 

S(y>,t)  =  w2(r+/zy>)  +Ri(t)w1(r+/xy>) 

(24b) 

-  ,  .  /xw'(r)+qw,(r)  d 

Ro(l)  =  -  ^(r)+qw2M  •  V2W=drw,2(r) 

(24  c) 

—  (‘wt(t-y|)w^(r+fiy,)+w;(t-yl)w2(r+^yi) 

*  ~  ^w1(t-yi)w1,(r+/iyi)+w'(t-yi)w1(r+/iyi) 

(24d) 

and 

t-yi  =  /i2(r+/iy,) 

(24e) 

Inserting  the  expansion 

(1-R)-1  =  SR"  .  R--=R0(t)Rj(t)  , 
n=o 

(25) 

into  (24),  and  Interchanging  the  order  of  summation  and  Integration,  one  may 
order  the  resulting  series  as  follows, 

4  00 

*  =  £  tf())  ,  *(i)  =  SGnj)  (26) 

i— i  n=o 

where  G^,  J=1...4,  represents  four  species  of  generalized  ray  Integrals  (Fig.  4) 
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Gn1}  =  ~f  ^ i ( r +/* y< ) w2( r -f a* y > )R Vxtdt 

4f*  c 

(27a) 

G£2)  =  TT  /  w2(r+/iy<)w2(r+/iy>)S0Rn5ixtdt 
c 

(27b) 

G^3)  =  -/  w1(r+Ai<)w1(r+^y>)RiRnelxt,dt 

4/i  c 

(27c) 

Gn4)  =  T~  I  w2(r+/iy<)wi(^+hy>)R0RiRneixtdt 

4/i  c 

(27d) 

The  physical  content  of  the  generalized  ray  fields  may  be  understood  by 
employing  asymptotic  approximations  in  the  integrands,  and  then  evaluating  the 
integrals  by  the  saddle  point  method.  For  large  values  of  their  arguments,  the 
Airy  functions  wt  2  in  (23)  can  be  approximated  by  [5] 

w1>2(-  u)  ~  T  br“1/V‘1/‘1  exp 

wl<2(t/)  *P  l7r"1/,2i/“I/4exp| 

Propagating  ray  fields  are  established  by  progressing  phase  terms  as  in  (28a). 
Because  of  the  different  arguments  of  the  various  w1  2  functions  in  (27),  one 
must  Identify  distinct  t-Intervals  where  these  arguments  have  positive  or 
negative  real  parts.  For  the  present  purposes,  to  obtain  propagating  ray  fields, 
it  is  required  that  a)  Re  (r+/xy<)>0,  b)  Rer<0;  but  one  may  have  c) 
Re(t-yj)<0.  As  seen  from  (24d),  the  argument  (t-yj  appears  in  the  reflection 
coefficient  due  fo  the  permittivity  slope  discontinuity  at  y=yt.  To  distinguish  the 
two  cases,  we  append  the  subscript  (1=1  for  Re(t-yj)>0  and  (1=2  for 
Re(t-yi)<0.  The  ray  integrals  pertaining  to  the  spectral  intervals 
defined  by  the  constraints  a),  b)  and  ^  2)  are  then  reduced  to  the  following 
forms  (see  Fig.  5  for  typical  trajectories): 


±  iJLj/3/2  -tiTr/4,  ,  Reu>0  (28a) 

3 

— t'3/2  ,  Re  r>0  (28  b) 

3 
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j=L_L  1=  1 

Here,  via  (28b),  Rj(t)~-i,  which  implies  that  the  upgoing  waves  in  this 
spectral  interval  are  evanescent  at  y=y|.  Virtual  reflection  takes  place  at  a  level 
y<yif  and  the  value  (-1)  of  the  reflection  coefficient  Identifies  an  encounter  with 
a  caustic.  Then, 


Gb=' (-l)"rn(t;p(0);q) 


exp{ix*|j!(t)} 


47 r 


{p(y<)-t}1/4{p(y>-t}1/4 


(29) 


where 


--sip 


(29  a) 


*(3(0  =  |2n(p(0)- 1)3/2  +  (p(y<)- 1)3/2-  (p(y>)-  t)3/2J  (29b) 


The  forms  for  J=  2,3,4  are  similar.  The  saddle  point  t$3  of  the  integrand  in 
(29)  is  located  from  the  stationary  phase  condition 


d_ 

dt 


*{3(t)  =  °  at  t[3 


which  yields  the  solution 


(30) 


•-  ~  (2n(p(0)-t)l/2+(p(y<)-t)1/2-(p(y>)-t)l/2|  (31) 

/i3  1  tfja 


Then  by  the  conventional  saddle  point  technique  [6], 


(-  i)nrn(t(3;p(o);q) 


exp  (ix<E  ( t{3 )  +( ivr  /4  )sgnQ|  ^  } 

{p(y<)-t|3}1/4{p(y>)“  ^S)1/4 

(32) 


with 


/i3X 


2n 


2Vp(0)-t  2v/p(y<)-t 


1 

2V/p(y>)~ 


(32a) 


The  saddle  point  condition  in  (31)  Is  the  equation  for  a  ray  path  reflected  n- 
tlraes  at  y=  0,  the  term  Tn  In  (32)  incorporates  the  associated  reflection 
coefficient  at  y=  0,  the  exponential  term  in  (32)  describes  the  phase 
accumulation  along  that  path,  (-l)n  accounts  for  n  ray  encounters  with  a 
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caustic,  and  the  remaining  amplitude  terms  account  for  ray  tube  spreading. 


Now,  from  (28a) 


R(t) 


if  (t;/x;yj)e 


3  na 


(33) 


where 


WgiJ/Wg^fr)} 

A»+MA*8€i)/wx(A*aW}  {wx(£j)/w '(£,)} 


(33  a) 


Then 


gO) 

u2,n 


/  .rn(t;p(0);q)fn(t;Ai;y|) 

47r  ca 


exp{lx<fr2(In)(t)} 
(p(y<)-  t}1/4{p(y>)- 1}1/4 


(34) 


where 

*&>(»)  =  w--^-  {2n  [(p(o)-t)3/2-(p(y,)-t)3/2]  +  (p(y<-t)3/2- (p(y>)- t.)3/2j 

(34a) 


The  saddle  point  condition  yields 

x  =  ~{2n[(p(0)-t)l/*-(p(y,)-t)1/4]  +  (p(y<)-t)1/2-(p(y>)-t)1/2J  (35) 
and  the  asymptotic  result 


1_ 

4 


•rn(t|,1IJ;p(0);q)fn(t|>1I);/i;yi) 


exp{lx»^1n)(t|1n))+(i7r/4)sgnQ^1n)} 

{p(y<)“  4.n}l/<{p(y>)-  ^.n)1  M 

(30) 


1 

lt3x 


1 

2\/p(0)-t 


2\/p(yi) 


1 _ ' 

(yO-t. 


+ 


A 


2  >/p(y<)-t  2  v/p(y>)-t  j 

(36a) 


Similar  results  can  be  derived  for  the  other  species  J=  2,3,4,  $=I,2. 

The  forms  of  the  solution  for  #=  1,2  and  J=l  in  (32)  and  (30)  are  generic 
also  for  the  other  ($tJ)  combinations.  Thus 
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rn(^;p(0);q)fn(^;/i 


;yi) 


exp[lx<l»^)(^)+(l7r/4)sgnQ(j^] 

(37) 


n  (y,t)  =  [p(y)~  t] 1/2 


(37  a) 


Here,  r  is  defined  in  (29a),  and  f  in  (33a),  for  5=2,  whereas  F=-i  for  5=1* 
Moreover, 


=  ^j>+^-{(2n+ai)  [a2n  3(0,t)-  a3n  3(yi,t)]+a2/n  3(0,t) 

-  a£n  3(yjft)+a4n  3(y<,t)+a5fi  3(y>,t)| 

>  tj!l 


(37b) 


Q®=  2^{(2n+ai)  t®*0 " 1(0,t)”  a®n  "  1(yi't)]+a2n  "  '(Oat)-  *£n  -  l(yittj 


+  a4n-1(y<»t)+as^_1(y>*t)J 


t# 


and  the  ray  equation  (saddle  point  condition)  is 


(37c) 


x=-^-J(2n+a1)  [a2n  (0,t)-  a3n  (y,,t)  j+a^n  (0,t)-  a^Q  (y^t)  (37d) 

4-  a4n  (y<,t)+asn  (y>,t)[  (37d) 

>  t U 

Then  the  following  definitions  of  the  constants  av  i=  l  to  5,  apply: 

{=1,  J=2:  a1—a2=-  a4=-  a5=l,  a,j=a^=a3=0  (38a) 

5=2,  J=2:  a i =<23=0,  a2=2,  a3=- a4=- o;5=l  (38b) 

5=1,  J=3:  a1=a3=o:2=a3=0,  a2=a4=as=l  (38c) 

5=2,  J=3:  oii=a2=0,  a2=o:4=a5=l,  a3=2  ^  (38d) 

5=1,  J=4:  a1=2,  a2=- a4=a5=l.  a3=a2=a^=0  (38e) 

5=2,  J=-l:  a1=2,  a2=a3=- a4=as=l,  a2=Q;3=0  (3Sf) 


Typical  ray  paths  are  identified  in  Fig.  5.  All  of  the  asymptotic  ray  field  results 
above  are  based  on  the  valaidity  of  Hu*  asymptotic  approximations  in  (2S)  for 
the  integrand,  and  also  on  the  validity  of  the  isolated  saddle  point  evaluation. 
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This  excludes  observation  points  near  cautlcs,  and  also  near  glancing  rays  that 
are  tangent  to  the  refractive  index  slope  discontinuity  (duct  boundary)  at  y=y, 
(Fig.  5).  _  ‘ 

D.  Hybrid  Ray-Mode  Expansion 

The  most  general  format  Involving  ray  fields  and  mode  fields  is  obtained  by 
combining  these  self-conslstently  in  a  hybrid  ray-mode  expansion  [4].  To 
generate  the  hybrid  representation,  we  first  deform  the  integration  path  C  in 
(19)  Into  the  two  contours  C^,  £=1,2,  which  surround,  respectively,  the  trapped 
modes  with  Retm>p(yj)  and  the  leaky  modes  with  Retm<p(yj);  these  contours 
(Fig.  8)  establish  the  rigorous  complex  extension  of  the  real  spectrum  intervals 
associated  with  £=1  and  £=2  In  Section  W.C.  Next,  we  expand  the  resonant 
denominator  not  into  an  infinite  (ray)  series  as  in  (25)  but  into  a  truncated 
series: 


l 

1-R 


-JL.  -  2.  i±* RN’  -  i-R*’  +  Rn  +  —  -^±4rN' 
1-R  2  i-R  2  n±^a  2  1—  R 


(39) 


Then  decomposing  the  attentuation  function  corresponding  to  and  C2  as  in 
(26),  one  obtains 

(40a) 

l-i  9=i 

where 


Here,  are  the  generalized  ray  Integrals  associated  with  Cp,  £=1,2,  which 
have  the  Integrands  shown  In  (27).  The  remainder  Integrals  R(J$j,  have  for  each 
wave  species  J=  1  to  4  the  same  integrand  as  the  corresponding  ray  integral 
Gg  N(,  but  contain  in  the  Integrand  the  following  additional  factors: 

R^i,:  M,  M=~(l+R)(l-R)"1  ,  J=1  to  4  (41a) 

2  s 

R0,,:  (RNa(l-R)]-1  -  M  ,  J=1  to  4  (4lb) 

Asymptotic  reduction  of  the  ray  Integrals  proceeds  as  in  Section  IV.C.  The 
remainder  integrals  can  be  reduced  by  deforming  the  integration  paths  C1<2  into 
steepest  descent  paths  SDP^  through  the  saddle  point  t£Jf  of  the  corresponding 
ray  integral.  Because  each  integrand  contains  the  resonant  denominator 
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(l-R)-1,  some  of  the  modal  pole  singularities  are  encountered  during  the 
deformation  (Fig.  7)  and  contribute  residue  contributions  that  furnish  partial 
modes  ("partial”  implies  that  the_  decomposed  traveling  wave  spectra 
represented  by  wx  and  w2  in  the  integrands  contribute  individually,  instead  of 
the  full  spectra  expressed  by  the  Ai,  Bi  functions  In  (22)).  For  $=J=1  and  Nt 
one  obtains  from  deformation  of  Cx, 

S  +  («) 

m=l 


where 


wi(rm+^y<)w2(rIn+^y>)  ix^ 

U. ' 


(42a) 


is  the  partial  mode  field,  and  r£M)(Ni  Is  a  new  remainder  integral  having  the 
same  integrand  as  R{$t  but  extending  along  SDPj$,.  For  saddle  point  tNi 
separated  from  the  modal  pole  t^n,  the  new  remainder  integral  can  be 
evaluated  by  the  saddle  point  method  to  yield 


(43) 


where 


L  = 


l+RQ(t)R|(t) 

l-R0(t)Rj(t) 


I  {■!<-  r)^»+x/4} 

e  3  -  e  3  T 


e  3  +  e  3  r 


(43a) 


and 


r  =  r(t;p(0);q)  =  - 


q-  l\/p(0)-  t 
q+f\/p(0)-t 


For  a  perfectly  conducting  earth  (q=oo),  (43)  reduces  to 


-7G[.Nt  cot  tX.1^,)3/2+7r /4) 


(43b) 


(44) 


The  contribution  in  (43)  or  (44)  is  usually  small.  When  the  saddle  point  is 
near  the  modal  pole,  the  new  remainder  integral  cannot  be  neglected  and  must 
be  evaluated  by  uniform  asymptotics  In  terms  of  incomplete  error  functions  [). 
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For  #=J=1,  and  N2,  the  remainder  integral  is  split  into  two  parts  (see 
(41b)).  For  the  first  term  in  (41b),  the  integration  path  Cx  is  deformed  around 
the  enclosed  pole  singularities  to  generate  a  sum  of  modes  as  in  (42a).  For  the 
second  term  in  (41b),  Ct  is  deformed  into  SDP|$a.  This  yields 

R{&  =  E  <4° -  R,  .  as Mj&  (45) 

m=q 

The  integral  Rq  is  the  same  as  the  one  for  r|$2  arising  from  the  second  term 
of  (41b),  except  that  the  integration  path  is  SDP[^2.  Therefore,  the  evaluation 
of  the  integral  can  proceed  as  above.  The  same  scenario  can  be  repeated  for  all 
of  the  other  integrals  R0j,.  Details  are  given  in  Appendix  A.  When  all  of  these 
results  are  combined,  one  obtains  the  following  most  general  hybrid  ray-mode 
representation  of  the  attenuation  function: 


V»=  E 


j=i 


4  «©,  4  M&,  4  CO 

+  E  S=?  +  S  2  4"  +  S  E  4P 

j=l  m=l  J=1  1=1 


+  E  E  PmJS,  -  RM$J  (46) 

i=i  5=i 

This  expansion  is  exact  if  the  generalized  ray  integrals  and  all  remainder 
integrals  are  kept  Intact,  but  useful  results  are  obtained  by  employing  the 
asymptotic  reductions  discussed  earlier. 
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Fig.  1  Spherical  earth  (radius  r=  a)  with  vertical  electric  dipole  source  at  Q. 

(r,0,^)  are  spherical  coordinates,  and  the  excited  field  components 
are  Er,  Eg,  H^. 


Fig.  2  Bilinear  permittivity  profile  corresponding  to  (10)  and  (14). 
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Fig.  3  Qualitative  placement  or  singularities  (for  low  loss  surface 
conditions)  and  integration  path  for  (19)  in  complex  t-plane. 
Trapped  and  leaky  mode  ranges  correspond  roughly  to  Retm>p(y,) 
and  Retm<p(yj),  respectively,  where  y4  locates  the  upper  boundary 
of  the  duct.  m=  mode  index. 
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\j=3  j=4 

'earth's  surface 


Pig.  4 


V/  1-4 1  VII  V/  WIV4I  I  VA  W  *-» 

(b)  1=2 

Ray  species  J=l  •  •  •  4  for  categories  $=l  (trapped  Inside  duct)  and 
$= 2  (leaking  out  of  duct).  For  arbitrary  reflection  index  n,  the  ray 
species  are  organized  according  to  their  directions  of  departure  from 
Q  and  arrival  at  P. 


leaky  ray  /  /  /  fy  t  duct  boundary 

\  ./x  /  J  // 

yi-  -  >7<  "  /l>^(Un)=(3,1,2) 


/ 


■(j,l,n)=  (3,2,3) 


/  trapped  rays 
glancing  ray 


earth's  surface 


Fig-  5  Typical  ray  trajectories  categorized  as  in  Fig.  4.  Trapped  rays 
leaky  rays  glancing  rays  — . 
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Fig.  6  Deformed  Integration  contours  preliminary  to  hybrid  ray-mode 
expansion.  and  C2  enclose  the  trapped  and  leaky  modes, 
respectively. 


Fig.  7 


Hybrid  ray-mode  expansion  by  deformation  of  integration  contours 
Cj  and  C2  into  steepest  descent  paths  SDPjj^  through  the  saddle 
POims  Ox  tis€  remainder  Held  integrands.  Mode  contributions 
arise  from  poles  crossed  during  the  deformation. 
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DOCTRINE. C  MAINTENANCE  DESIGN  NOTES 

The  AEGIS  doctrine  parser  is  designed  to  provide  a  high 
degree  of  operational  flexibility  without  burdening  the  human 
operator.  This  flexibility  could  easily  permit  a  DOCTRINE 
derivative  to  be  used  to  extract  information  of  any  potential  use 
as  long  as  the  information  is  preceded  by  a  well-defined  trigger 
phrase. 

The  syntax  to  use  DOCTRINE  follows: 

DOCTRINE '[doctrine-file-name  [message-source  [trigger-phrase]]] 

DOCTRINE  is  intended  to  be  run  in  background  mode;  it 
requires  no  additional  action  on  the  part  of  the  user.  During 
operation,  the  program  scans  messages  arriving  over  a  serial  port 
until  AEGIS  doctrinal  information  is  received.  The  doctrinal 
information  is  then  stored  into  a  file.  This  process  continues 
indefinitely. 

To  make  operation  as  easy  as  possible  on  the  user,  default 
values  for  the  command  line  parameters  are  assumed  if  none  are 
provided.  The  following  defaults  are  hard-wired  in  function 
READ  ARGUMENTS  of  DOCTRINE. C: 

"destination”  —  Where  is  the  extracted  doctrine  information 
to  be  stored?  If  no  "doctrine-file-name"  is 
provided  on  the  command  line,  DOCTRINE  sends 
its  extractions  to  MM  (empty  string) . 

"msgsource"  —  Where  is  the  information  to  be  received?  If 
no  "message-source"  is  provided  on  the 
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command  line,  DOCTRINE  attempts  to  read  the 
"/dev/ttyb”  serial  port. 

’•trigger”  —  What  phrase  is  used  to  indicate  the  start  of 
AEGIS  doctrine  in  the  message?  If  no 
"trigger-phrase"  is  provided  on  the  command 
line,  DOCTRINE  uses  "XXXAEGISDOCXXX” . 

In  order  to  provide  the  greatest  flexibility,  DOCTRINE  » 
searches  for  occurrences  of  "trigger-phrase"  in  the  incoming 
message  stream.  This  not  only  avoids  the  slower  finite-state 
approach,  but  also  permits  any  number  of  ways  that  the  doctrine 
message  could  be  formatted  (’’NOTE”,  "NARRATIVE”,  "REPLY”,  "TEXT”, 
etc.).  This  disadvantage  to  this  approach  is  that  ANY  occurrence 
of  the  trigger-phrase  initiates  extraction.  Consequently, 

«  i  * 

messages  that  discuss  the  trigger  phrase  but  do  not  actually 
contain  doctrinal  information  should  not  be  passed  to  the  parser. 

The  particular  algorithm  that  DOCTRINE  uses  to  detect  the 
trigger  is  fast,  but  not  particularly  robust.  It  can  miss  some 
occurrences  of  the  trigger  when  the  trigger  is  preceded  by  a 
partial  phrase  that  looks  like  the  trigger.  Fortunately,  this  is 
not  possible  with  the  default  trigger.  However,  if  the  trigger 
were  "ABABAC”,  then  any  doctrine  information  following  the  line 
that  contains  "ABABABAC”  would  be  missed  —  even  though  the  % 
trigger  phrase  is  present.  The  scanning  algorithm  does  have  some 
robustness  so  that  inadvertently  duplicated  characters  would  be 

s 

ignored.  ( "XXXAAEEGGHSSDDOOCCXXX"  would  trigger  extraction  in 

the  default  case.) 

In  its  delivered  form,  DOCTRINE  discards  all  characters  that 
follow  the  trigger  phrase  on  the  same  line.  That  is,  the 
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doctrinal  information  is  presumed  to  start  on  the  line  FOLLOWING 
the  trigger  phrase.  If  the  information  that  immediately  follows 
the  trigger  needs  to  be  extracted,  that  can  be  accomplished  by 
commenting  out  the  logic  branch  in  function  READ_TTY  that  begins 
•’if  (trigstat==0) " .  (The  integer  "trigstat"  provides  an 
indication  of  how  much  of  the  trigger  phrase  has  been  read.  It 
is  negative  if  not  all  the  phrase  has  been  found.  It  is 
incremented  as  each  received  character  is  matched  against  the 
trigger  phrase  in  memory.  When  "trigstat"  becomes  0,  a  match  of 
the  trigger  phrase  has  been  found.) 


After  the  end-of-line  that  follows  the  trigger  phrase  has 
been  received,  DOCTRINE  stores  all  new  characters  received  into 
the  doctrine-output-file  until  a  slash  ("/")  character  is 
received.  The  assumption  used  is  that  no  slash  character  will  be 
embedded  within  the  raw  doctrine  information.  (If  this  turns  out 
not  to  be  the  case,  the  process  for  determining  the  end-of- 
doctrine  will  be  more  complicated.  For  instance,  a  particular 


character,  even  a  slash,  can  be  used  to  prefix  all  doctrine  lines 
in  the  formatted  message.  A  more-straight-forward,  but  less 
flexible,  approach  is  to  merely  copy  a  specific  number  of 
characters.)  With  the  slash  character  as  end-of-doctrine,  the 
output  doctrine  file  could  contain  a  formatting  keyword  (such  as 
"ENDAT”) .  The  end-of-doctrine  character  is  hardwired  into  the 
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By  adjusting  the  output  file,  the  input  source,  the  trigger 
phrase,  the  start-of-data  (end-of-line  by  default)  and  the  end- 
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of-data  (m/m  by  default)  definitions,  DOCTRINE  can  be  modified  to 
extract  any  type  or  amount  of  data  from  formatted  messages. 
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Approved  for  public  release;  distribution  Is  unlimited. 

The  views  and  conclusions  contained  In  this  report 
are  those  of  the  contractors  and  should  not  be 
Interpreted  as  representing  the  official  policies, 
either  expressed  or  implied,  of  the  Naval  Ocean 
Systems  Center  or  the  U.S.  Government. 


